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ABSTRACT 

We study the production of {t + t)g at the hadron colliders in an R-parity 
{Rp) violating supersymmetric model. This process provides us with infor- 
mation not only about Rp violation, but may also help us in detecting the 
supersymmetry itself. It is possible to detect an Rp violating signal (with 
single gluino production) at the future hadron colliders, such as Fermilab 
Tevatron Run II or CERN Large Hadron Collider (LHC), if the parameters 
in the supersymmetric ]}ip interactions are not too small, e.g. for rrig = 1 
TeV, a" =0.1, still hundreds of events are produced at LHC with luminosity 
30 fb~^. Even if we could not detect a signal of ^p in the experiment, we 
get stringent constraints on the heavy flavour ^p couplings. In addition to 
the minimal supersymmetric standard model we have also considered some 
models with a heavy gluino as the lightest supersymmetric particle. 
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I. Introduction 



It is well known that the minimal supersymmetric standard model (MSSM) Q in its most 
general form contains lepton number (L) and baryon number (B) violating couplings. The 
resulting catastrophic proton decay Q can be avoided by imposing R-parity symmetry Rp, 

= (_l)3B+L+2S^ (1) 

where S is the spin of the particle. In the models with R-parity conservation, superparticles can 
only be pair produced and the lightest supersymmetric particle (LSP) will be stable. However, 
Rp conservation is not necessary for forbidding proton decays, instead of that, we just need 
either B-conservation or L-conservation Models with Rp violation (^p) can provide many 
interesting phenomena, such as neutrino masses and mixing. Partly because of that, ]^p has 
attracted much attention Q |5|, |6|. Many constraints from low-energy phenomenology are 
cohected in Ref. [0]. 

The supersymmetry (SUSY) must be broken, since it has not been observed so far. Two kinds 
of breaking mechanisms of supersymmetry have been extensively studied phenomenologically, 
namely the minimal supergravity (mSUGRA) and minimal gauge-mediated SUSY breaking 
(GMSB) Q. These predict a different pattern of masses especially for the partners of gauge 
bosons, gauginos. Thus, finding the signals of gauginos is an important way to probe SUSY. In 
mSUGRA and GMSB, it is assumed that masses of gauginos will be unified at Grand Unified 
Theory (GUT) scale. From the evolution of parameters, gluinos should be the heaviest gauginos 
at low scale, since the ratios of gaugino masses to coupling constants do not change with scale 
in one-loop approximation |jlO[. Thus gluinos can decay to other gauginos with jets in the 



i?p-conserving model [11|. 



However, heavy gluino as LSP (or NLSP) may still exist in some GMSB models, as S. Raby 
suggested in fl^ . He introduced the Higgs- Messenger mixing in GMSB model. This would lead 
triplet-doublet messengers split from triplet-doublet Higgs splitting. Since gluino obtains mass 
from SUSY-breaking induced by triplet messengers, heavier triplet messengers could suppress 
the mass of gluino so that gluino could be the LSP (or NLSP). The experimental limits on 



the masses of gluinos have been discussed in i?p-conserving model |13]. In i?p- violating model 



gluinos can decay through i?p-violating channels, which obviously changes both the detection 
strategy and present mass limits [p^ . 

In the high energy hadron colliders especially squarks and gluinos will be produced plentifully, 
if they are not too heavy. It is hoped for that information on SUSY is found in Fermilab Tevatron 
Run II, or in the future CERN Large Hadron Collider, where squarks and gluinos with masses 
below 1.5 TeV should be detected In mSUGRA and GMSB models, masses of squarks and 
gluinos are usually of the same order. However, in some special mechanism, such as O-II model 



|16|, it has been suggested that squarks can be much heavier than gauginos. In this model the 
gluinos will be produced at much lower energies than squarks, so that it can be the first detected 
SUSY particle at hadron colliders. Since gluinos would be almost degenerate with the lightest 
neutralino and chargino, Rp violating decay of gluinos will be significant 0. 

The single production of gluinos, neutralinos and charginos has already been considered in 
the general case [p!7|], and the single production of squarks, which is also significant to detect 
SUSY and i?p- violation, has been considered recently In this work we will consider ^p at 



hadron colliders in the process 
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This process occurs via B-violating terms in the Ifip model. In terms inducing heavy flavours, 
the-^p couplings can be very large from the present upper limits 0]. For example, \2ij and X-nj-, 
getting their strongest constraints from the ratio of widths of Z to leptons and hadrons, can be 
of order one (0(1)) for the sfermion mass 0(100 GeV). 

The pair production and decay of gluinos at hadron colliders have already been researched 
m Ref. It has been shown that detecting gluinos is very difficult. It may become easier 

with an accompanying top quark in the process which we consider. 



On the other hand, the single top quark production is an interesting topic itself |2C]. The 



single top quark production in the I^p model has been considered in Ref. [21|. There, the 
possible cross section will depend on the I/ip parameters as |A |^, while the process we will 
consider depends on |A"p. Although the mass of the gluino is not known presently, it is still 
possible to get stronger constraints on iZp-violating parameters from the {t + t)g production. 

In the following, we will give the analytical calculations of PP{PP) tg{tg) in section II. 
In section III gluino and top decays are considered and in section IV the numerical results are 
presented. The conclusions are given in section V and some details of the expressions are listed 
in the appendix. 

II. Production of gluinos in PP — > tg 

The superpotential for violating, but gauge and supersymmetry preserving interactions 
is written as 

^i/ip = \ij]kU-LjEk + X[jkLi-QjDk + X"[ji,fiiDjDk + eiLiHu (2.1) 

where Lj, Qi and are SU(2) doublets containing lepton, quark and Higgs superfields respec- 
tively, Ej {Dj, Uj) are the singlet lepton superfields (down-quark and up-quark). The square 
brackets around the generation indices i,j denote antisymmetry of the bracketted indices. 

We ignored the last term in Eq. (2.1) because its effects are assumed small in our process ptl. 
Thus, we have 9 A-type, 27 A -type and 9 A -type independent parameters left. The constraints 
on the couplings 



m \ 
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|(AorA)A I < 10-^*^ ( . (2.2) 



100 



is usually taken to indicate that only L- or B-number violating couplings exist. 

In our work we will only consider the baryon number violating couplings, i.e., the third term 
in Eq. (2.1). In the following calculations we assume the parameters A to be real. 

We define the Mandelstam variables as usual 

s = ipi+P2f = {P3+P4?, (2.3.a) 

t = {pi-P3f = {Pi-P2f, (2.3.5) 

u = {pi- Pif = {P3 - P2f. (2.3.c) 
The amplitude (Feynmann diagrams in Fig.l) of qjq'k tg is given by: 

M = Ms + Mt + Mu, (2.4) 
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+ crossed diagram 



Figure 1: Feynman diagrams of qjqk tg. 

with 

Ms = lF'{ps){-iV2gsT^^> ) (cos OPl + sin ePR)va' (Pi) 



+ u-ips){-iV2gsT^^,)ism9PL-cosePR)v^,{p,) ^^_^, ^'-^ 
X U^{p2){-2ie"^Px;,.)cos9PRUp{pi), 

Mt = i^(p3)(iV2ff.r^"^Pij)t.p(pi)(^<(p4)(-2ie"^^A;;.,Pfi)n^(p2), (2.5.5) 
Mu = ^(p3)(iV25sT^"^Pii)u^(P2)(^^<(p4)(-2ie"^^A3,^.Pfi)^Xp(pi), (2.5.c) 

where Pl,r ai"e left- and right-hehcity projections respectively, is the mixing angle of stop 
quarks (see Appendix for details) and F^^ ^ ^^e decay widths of stop quarks ti^2- The families of 
down-type quarks are marked by j, /c = 1, 2, 3 and the upper index c means charge conjugate. In 
the calculations, we have neglected all mixing angles of scalar quarks except stop quarks. The 
amplitude depends on the i?p-violating parameters X^^^^ {j, k = 1,2, 3) in the process. 

The total cross section for the process qjq'f^ tg is: 

<t(s) = j^j^ diEspins \M\^ (2.6) 



where = ^ 



[ml + m| — s) ± w + + m| — 2smf — 2sm| — 2m|?n,| 



and M is the ampli- 
tude. Here we have neglected the masses of incoming quarks. A'^c = 3 is the color factor and the 
bar over summation means averaging over the initial spins. 

In a similar way, the cross section for g^^^ — > tg can be calculated. The possible effects of 
^3% ~^ ^9 ^^'^ ~^ ^9 should be observed in PP or PP colliders. The cross section for the 
process P{Pi)P{P2) — > (t + t)gX can be obtained by convoluting the subprocess with quark 
distribution functions [p^ ], 

a{s) = J dxidx2fi{xi,Q)fj{x2,Q)a{s,asifJ.)) (2.7) 

with pi = xiPi, T = xiX2 = s/s. fi.j{xn, Q) (n = 1, 2) are the corresponding quark distribution 
functions of protons. We take Q = fi = 300 GeV. Similarly we can find numerically the cross 
section of P{Pi)P{P2) ^ {t + t)gX'. 
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Figure 2: ratio of Fj^^ : Thiy as function of mass of stop quark ti, with nig = 30 GeV. 

III. Top and gluino decays 

In the process of Eq. (2), we have in the final state two heavy particles which will possibly 
decay inside the detector. Here we will shortly review the relevant decay modes of the top quark 



and gluino p|, p|, p|, KM 



The top decays in the MSSM have been considered by several authors ||2^, |2j, |25|, P?]] . 
The main decay mode is t — > bW, but t bH~^ can compete with it if mass of the charged 
Higgs is lighter than mt — m{,. Top quark decay to R-odd particles will also be important if 
those superparticles are light enough |2^. However, in our case, with heavy squarks, the decays 
to real superparticles are impossible, except for the light LSP gluino with light squark [^]. We 
plot the ratio F^^ - : F{,vk as a function of the stop quark ii mass, m^^ , in Fig.2 for the gluino mass 
rrig = 30 GeV. In order to guarantee the purported standard top quark events at the Tevatron, 



BR{t — > hW) should be larger than 40 — 50% as lower bound [28|. So with the assumption of 
light LSP gluino (about 25 — 35 GeV), lower limit on the mass of stop quark can be obtained. 
For the top quark decay through iip- violating interactions [^5|, the branching ratio of those 
decay modes will be very small compared with t — > bW in our case. The only decay channel 
which we will consider in detecting top quark is the t bW — > blvi, where I = e,^. We will 
confine us to these decays (with BR(W lui) ~ 22%, its branching ratio in top quark decay 
should be at least 8.8 — 11%) since it is assumed that they have less background and are thus 
easier to detect than the hadronic decay modes. 

We consider in this work a heavy gluino, which may be the LSP. The decay modes of a heavy 
gluino have been looked at e.g. in |14, 15| in Rp violating case. The decay channel g ^ qq will 



dominate if kinematically allowed. A lighter than squarks gluino will decay to 

9 ^ qqXi, m'xf, 9Xi, (3) 

where and xf are neutralinos and charginos, respectively. The R-parity breaking decay 
modes become important for large A": 

9 qiqjqk- (4) 

In Fig. 3, we draw the g decay with parameters in mSUGRA model. We consider the decay of 
gluinos in Fig. 3 (a) and (b) with Rp conservation and violation, respectively. In our calcula- 
tions, only the two lightest neutralinos and the lightest chargino are considered, since the other 
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Figure 3: a) Decay width of gluino as a function of mass of gluino, rrig, for Rp conservation. 
Solid line corresponds to the total decay width and dashed line to partial width containing heavy 
quarks (top quark and bottom quark), b) Decay width of gluino as a function of rrig for Rp 
violation. Solid line corresponds to the total decay width and dashed line to the partial width 
containing heavy quarks (top and bottom), c) Decay width of gluino to tti as a function of rrig. 

neutralinos and the heavier chargino are too heavy for gluino to decay into. It is shown that 
the gluino branching ratio to two heavy quarks (top quark or bottom quark) and neutralinos or 
charginos (or jets from i?p-violating interactions) can be very large and increases with mass of 
gluino. This is reasonable because stop quarks and sbottom quarks can be much lighter than 
the other scalar quarks due to the large Yukawa couplings. Especially, in i2p-violating terms, the 
process through virtual stop quark ti will dominate the decay width because in the mSUGRA 
model it is the lightest squark and its mixing angle is near tt/2. Thus, there are two b (or t) 
quarks as a signal of gluino in our process (branching ratio can be read from Fig. 3, as about 
60 — 70%). Combined with other top quark, the three heavy quarks, leads to the final state 
36 + n(Z + I') {n = 1, 2, 3). This final state can be detected in the future CERN LHC and dis- 
tinguished from background (with an assumed b-tagging efficiency e ~ 50% in LHC [§^]). For 
the much heavier gluinos which can decay directly to top quark and stop quark, decay width is 
shown in Fig. 3 (c). 

When gluino is the LSP, the only available decay modes are the R-parity violating ones. If 
the R-parity violating couplings are not exceedingly small, the gluino will decay in the detector 
through these channels. However, if the couplings are very small, the gluinos will form so-called 
R-hadrons before decaying (see e.g. [^, [T^). 

IV. Numerical results 

In the mSUGRA model we are interested in the region of the parameter space where the 
gluino is lighter than the squarks, with the possible exception of the lighter stop quark. This 
choice allows gluino to be produced via the third generation iZp-breaking coupling with relatively 
large cross section and on the other hand it is not complicated by gluino decay to other squarks 
than possibly the lightest stop. 

As a representative example of this part of the parameter space, we take niQ = 1000 GeV, = 
— 1000 GeV, tan /? = 10 and sign{ii) = +1. Varying mi/2 suitably gives us the relevant gluino 
masses. The resulting masses for supersymmetric particles with the varied mi/2 are listed in 
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Figure 4: a) Cross section of PP {t + t)gX as a function of nig with Ecms = 2 TeV. Solid 
line corresponds to Xy,^^ = 1.0 and dashed line to X^^^ = 0.1. b) Cross section of PP {t + t)gX 
as a function of nig with Ecus = 14 TeV, X!!^^^ = 0.1. 

Table iQ. 

Table. 1 MSSM parameters with mo = 1 TeV, Aq = —1 TeV, tan/3 = 10 and sign(fi) = +, 
units of values are GeV in the table. 



mi 


m-g 


mi. 




"^6i 


rriq 




m~± ~ m-o 


120 


308 


457 


813 


790 


1020-1040 


49 


95 


140 


360 


468 


828 


803 


1028-1055 


58 


111 


180 


463 


499 


862 


836 


1060-1090 


75 


144 


200 


515 


517 


882 


855 


1075-1110 


83 


161 


220 


566 


537 


904 


876 


1100-1130 


92 


177 


250 


643 


571 


939 


910 


1120-1160 


104 


202 


280 


721 


608 


977 


947 


1160-1200 


117 


227 


300 


773 


634 


1004 


974 


1190-1230 


125 


243 


330 


850 


675 


1047 


1016 


1220-1270 


138 


268 


340 


876 


689 


1061 


1031 


1240-1280 


142 


277 


350 


902 


703 


1076 


1045 


1250-1300 


146 


285 


370 


954 


733 


1106 


1076 


1280-1330 


155 


302 


390 


1006 


763 


1138 


1107 


1310-1360 


163 


318 



In Fig. 4 (a), we show the cross section of PP {t+t)gX as a function of mass of gluino (nig) 

at Tevatron Run II energy, i.e. with center-of-mass energy of the collision ^/s = 2 TeV. There, 

I'll 1 \ \ " \ " 11111 \" 

the solid Ime corresponds to X^i2 = A313 = A323 = 1 and dashed Ime to A312 = -^313 = -^323 = 0.1. 

At LHC with = 14 TeV, the cross section for PP {t + t)gX as a function of nig is shown 

in Fig. 4 (b) with X^i2 = X^i^^ = A323 = 0.1 for solid line and A312 = •^313 = -^323 = 0-01 for 

^We thank A. Wodecki for providing us his program, which calculates the sparticle masses in MSSM. The 
program checks also some phenomenological constraints. 
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Figure 5: Decay width of stop quark t2 as function of mass of t2- Solid line: X^^j = 1.0. Dashed 
line: X^^ =0.1. 



dashed line. In the calculations, we have neglected the decay width of ti, since mass of ti is far 
from the center-of-mass energy. T^^ including i2p-violating contribution is shown in Fig. 5, where 

solid line corresponding to X^i2 = X^i^ = A323 = 1 and dashed line to A3^2 — -^313 — -^323 = 0.1. 
The results show that the cross sections can be very large if nig is smaller than 400 GeV and A 
are close to the present limits (A" ~ 1) in the Tevatron Run II (with luminosity about 2fb~^, it 
corresponds to hundreds of events). At LHC with luminosity 30 fb^^, the process can potentially 
be seen with mass of gluino less than 1 TeV even if A < 0.01. 

We have also considered production of the LSP heavy gluino. It has been shown in the 
model of |13| that masses of gluinos below 115 GeV are excluded except a narrow window 



between mg = 25 — 35 GeV. We found that the cross section of the process in Tevatron Run 
I can be very large with masses of gluinos staying in that narrow window (about 10 pb with 
rUg = 30GeV,A3jj = 1 and all = 150 GeV). Unlike in the discussion of [^], gluinos can 

decay through i?p-violating interactions (e.g. g bcs) if A" is nonzero. The gluino decay added 
by single production of top-quark, should have been visible already in Tevatron Run I using 
above results. Therefore, the narrow LSP gluino window can give much stronger constraints on 
i?p-violating parameters A (with A =0.1, still several events are produced with luminosity 
19 pb~^ in Fermilab), otherwise we can close this narrow LSP heavy gluino window. 

In the mSUGRA model, the gaugino soft-SUSY-breaking masses are assumed to be equal 
at GUT scale and will lead to m?, : m-± : m^-o 7 : 2 : 1 at low energy. However, there are 



other mechanisms such as O-II model ]16[] in which ratios of masses can be given as 7713 : m2 : 
mi ~ —(3 + 6gs) '■ (1 — ^Gs) '■ (33/5 — 6gs) at GUT scale, where 5gs is the Green-Schwarz 
mixing term. Then at low energy, we can obtain light gluinos, which are almost degenerate 
with the lightest neutralino and chargino, with the heavy scalar quarks. As a representative 
of this model we take 6gs = —4.1, mo = 1800 GeV, mi/2 = 10 GeV, ^0 = 0,tan/3 = 3 and 
sign{n) = +1 (leading to rrig ~ 102 GeV and m^ ~ 1.5 TeV). Setting Ag^^ = 1, the cross 
section of the process PP — > (t + t)gX can be about 3.6 fb with y/s = 2 TeV and of the process 
PP {t + t)gX about 27 pb with -y/s = 14 TeV. So in this model the single production of 
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gluino can provide a very significant signal for detecting SUSY and Rp violation. Detecting pair 
production of gluino in this model was already discussed in |ll], |l5| and it was shown that the 
i?p-violating decay of gluino will dominate if A are close to the present upper limits. Thus, the 
production of heavy quarks from gluino decay may provide a good signal for detecting gluinos. 
Similarly, when we detect the production of (t + t)g, we have at least two heavy quarks for 
tagging. 

IV. Conclusion 

We have studied the processes PP[PP) — > {t + i)gX in supersymmetric models with explicit 
-Rp- violation. We have seen that it is possible to test the models at future Fermilab Tevatron 
Run-II and CERN LHC experiments, provided the A"-type couplings are large enough. We 
suggest also to check the old data in Tevatron Run I and get stronger constraints on Rp violation 
or exclude the narrow window of LSP heavy gluinos. 

The process, with single production of top quark and gluino, can give signals both for SUSY 
and Rp violation. Specifically, in a typical O-II model, the process will be an important one to 
check SUSY and Rp violation. It is also shown that gluinos should be detected through their 
i?p-violating decay if Rp violating parameters A are close to the upper bounds. 
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Appendix 

In this appendix we present the mass matrices jlO| j which we need in the calculations. 
A. The sfermion sector 

In the mSUGRA model, the masses of squarks and sleptons follow from the GUT scale 
parameters mo, ^o, tan /3 and sign{fj.) ||lO[. The sfermion mass matrices are given by 



M 



m 



|r + rnj mf{Af — firf) 



f y mf{Af — firf) m- +mj- 



(Al) 



where "T-j^ ^ are given in |1C], nif are the masses of the partner fermions and r^^s 



r, 



e,fi,T 



^/fu,c,t = tan/?. From these matrices we can get the mixing angles Of and masses of the 
sfermions: 

{A.2) 
{A.3) 



sm20f = — ' i " , cos26f = '-^ 

J m - —171. ' 

tl f2 



~!5 

m - —m- 

fl 12 



m 



fl,2 



m) + \ 



B. The chargino/neutralino sector 

In order to calculate the decay of gluino, we need to consider the chargino and neutralino 



mass terms. The general chargino mass matrix is given as follows |1C]: 

^ ( m2 V^ruu, sin /3 \ 
\/2muj cos /3 /i I 

It can be diagonalized by two real matrices U and V, 

U*McV-^, 

where U = 0-,V = 0+if detMc > and F = asO^ if detMc < 0, with 



0-3 



and 



+1 
- 1 

tan 26- 
tan 26^ 



cos 6± sin 0-1- 
- sin 6± cos 6± 



2\/2m™ (m2 cos f3+fi sin /3) 
m,| — /i^ — 2m J, cos /3 ' 

2-\/2»nm (»Ti2 sin /3+fi cos /3) 
mj— /i^+2m^ cos /3 ' 



Then the chargino masses are, 

x+ V2^ 2 I r I 



=F{("l2 



,2\2 



(-2 - H")" + 4m^ cos"' 2/? + Ami{m2 + ^ + 2m2//sin2/?)}2] 
In the case of the neutralinos, the mass matrix is given by |lO| : 

/ mi — rrizSu, cos P m^s^sin/? \ 

Mn = 



(B.l) 
iB.2) 

{B.3) 
(BA) 

{B.5) 



m2 rrizSw cos 13 —rUzSw sin (5 

—mzSwSm.fj nizCwSmfd — /i 

V rUzS^ sin (3 — rrizCwSinf] — fi 

It can be diagonalized by a unitary matrix (4 x 4) as: 



iB.6) 



{B.l) 
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